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PREFACE 


The objective of this investigation is to locate ocean current 
boundaries by sensing the color change associated with the cyclonic edge 
of the zone of maximum horizontal velocity shear. The test site is the 
eastern Gulf of Mexico where the strongly baroclinic flow from the Yuca- 
tan Straits forms into the Loop Current. The research is using ERTS' 
data in the investigation of ocean color sensing. 

A time series of the. location of the Loop Current was begun in 
August 1972. For every other transit of ERTS (i.e., 36 days) the current 
was located by ship. These data have been used to provide surface mea- 
surements in support of the spacecraft observations and to obtain the 
first time history of the circulation in the eastern Gulf of Mexico. 

These research cruises were conducted for over one year in order to 
understand the seasonal variability associated with the circulation and 
detection of the circulation's surface indicators. 

Surface chlorophyll-a concentrations, which contribute to the shift 
s in color from blue to green in the open sea, were found to have high 
spatial variability; significantly lower concentration were observed in 
the Gulf Loop Current. The cyclonic edge of the current is an accumula- 
tion zone which causes a peak in chlorophyll concentration. The dynamics 
of these currents also cause surface concentrations of algae, which have 
a high reflectance in the neaip infrared. Combining these observations 
gives rise to an "edge effect" which can show up as a brig.ht lineation 
on mul tispectral imagery delimiting the current's boundary under certain 
environmental conditions. Frequently, the sea-state in the current is 

i 
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higher than in surrounding water due to differential current shear. 

When high seas introduce bubbles, white caps, glitter, and foam, the 
reflectance is dominated by scattering rather than absorption. This has 
been detected in the imagery and used for current location. 

The color effect relates to the spectral variations in the optical 
properties of the water and its suspended particles, and is studied by 
radiative transfer theory. Significant oceanic parameters identified 
are: the probability of forward scattering, and the ratio of scattering 

to total attenuation. Several spectra of upwelling diffuse light are 
computed as a function of the concentration of particles and Gelbstoff. 
These calculations compare favorably with experimental measurements and 
show that the ratio-of-channels method to determine ocean properties 
gives ambiguous interpretative results. 

Subtle contrasts and low radiances observed by the ERTS multi spec- 
tral scanner over the ocean require computer enhancement for adequate 
analysis. Experiments designed to evaluate contrast stretching, ratio- 
ing, differencing, smoothing, filtering, and false color enhancing, in- 
dicate that the best information can be extracted by simple contrast 
stretching. Spectral analysis of the data shows that a low pass, two- 
dimensional filter kernel, designed to be 3 db down at 10 scanspots, 
effectively eliminates the six-line banding caused by the mul tispectral 
scanner design. Automatic contouring techniques require careful scrutiny 
because data fields are created which can lead to false interpretations. 
Joint histograms of oceanic radiances did not prove to be useful due to 
the low range of energy in the several spectral intervals. Comparisons 



of satellite data with surface ship observations confirm theoretical 
predictions of the difficulty in interpreting scenes of the coastal zone. 
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(ERTS ID 1026-15230). The Florida Current can be 
seen as a line of dark lineation parallel to the 
coast; bottom is essentially invisible in this MSS 6 
scene. A ship can be seen by its characteristic V- 
shaped wake just offshore of Miami Beach. Possibly 
the Virginia Key sewer outfall area can be observed 
by its low reflectance due to an organic slick. 

Horizontal distance across the image is 90 kilometers. 44 

Solid line: Along scanline power spectra of oceanic 

data in MSS 5, Dashed line: Cross scanline power 

spectra of oceanic data. Power is in units of DN 

which has a radiance equivalent of' mW cm" 2 sr _1 . 48 

d ... : . 
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Percent transmittance for the apex of the New York 
Bight. Upper [a]: 29 August 1973 (ERTS ID 1402- 
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in 8 [d]. Lower [b] : 16 September 1973 (ERTS ID 

1420-15067) transit processed as in [a]. Oceano- 
graphic stations used are located by dots. Nega- 
tive images. Horizontal scale in each figure 
approximately 50 kilometers. 


HISTORICAL INTRODUCTION 


Early work in the Gulf of Mexico has been summarized by Galtsoff 
(1954) who edited a complete overview of. the biology, chemistry, geology, 
and physics of the area. In 1895, Lindenkohl published a map of the 
temperature field at 250 fathoms. (457 meters), in degrees farenheit, 
which is reproduced in figure 1. The warm waters of the Caribbean can 
be seen flowing northward into the Gulf and penetrating deeply into the 
ambient thermal field; similarly the Gulf Stream, seen as a region of 
large horizontal temperature gradient, is flowing easterly and then 
northerly through the Strait of Florida. In an analysis of these data 
Sweitzer, in 1898, reported that the circulation was a spreading of this 
inflow which resulted in an anticyclonic circulation around the entire 
Gulf basin. Parr, in 1935, reported the opposite conclusion using 
ATLANTIS data taken in 1933: the Gulf Stream takes only the shortest 

path from Yucatan to the Strait of Florida. Even after reviewing 
Dietrich's 1939 map of the salinity maximum core,. which reflects the 
ydeep penetration seen in figure 1, Leipper expressed this divergence of 
opinion as the state of knowledge in 1954. Work done in the 1960's, 
notably by Leipper, Cochrane, Nowlin, and others at Texas A&M University, 
led Leipper (1970) to speculate that there was an annual cycle in the 
current patterns in the eastern Gulf, A time-series to establish the 
seasonality was reported by Maul (1974). This paper is a preliminary 
descriptive report on the time history of the current from August 1972 
to September 1973. 




Figure 1. Lindenkohl's 1895 map of the temperature field (°F) at 250 

fathoms in the Gulf of Mexico. Data are from soundings made 
by Sigsbee, Bartlett, and Brownson between 1874 and 1883. 

The shallow waters (<250 fms) of Campeche Bank and the west 
Florida Shelf, strippled on the chart, outline the topo- 
graphic constraints of the Gulf Loop Current. 
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OBSERVATION PROGRAM 


The field work began in June 1972 with a ship/aircraft experiment 
designed to detect the color boundary of the Loop Current front south of 
Dry Tortugas. The NASA C-130 flew over the research vessel track which 
was oriented in the same azimuth as ERTS suborbital track. Seven al- 
titudes were flown, at 100 mb decrements, over the ship. Aircraft data 
collection included RC-8 color and color IR photography as well as multi- 
spectral photography to simulate the ERTS MSS imagery, PRT-5 sea surface 
temperature profiles, Bendix 24 channel scanner data, and inflight re- 
cordings of atmospheric temperature, pressure and moisture. 

Prior to the aircraft overflights the scientific crew aboard ship 
located the boundary. Measurements of ocean temperature, chlorophyll -a, 
volume scattering function, and salinity were made every ten meters down 
to 50 meters; Fore! color, and upwelling and downwelling spectral irra- 
diance were taken at each of the five stations that bracketed the front. 
During the overflights on the following day, closely spaced surface 
measurements of the same variables (except irradiance) were taken. All 
measurements were made using standard oceanographic techniques or anal- 
yzed by methods detailed by Strickland and Parsons (1968). 

The sequence^ of pathlines of the Gulf Loop Current was started to 
obtain the location of the cyclonic edge by ship in order to evaluate 

several methods of location by satellite. Figure 2 is a cross-stream 

d . 

transect taken on May 1, 1973, from north to south along the suborbital 
track of the Earth Resources Technology Satellite (ERTS). Speed of the 
vessel was about 10 knots, so the horizontal distance traveled from 1200 

- 3 - 


y 


ORIGINAL PAGE 
OF POOR QUALU* 




Figure 2. Thermal cross-section of the Gulf Loop Current on 1 May 1973 
northwest of Dry Tortugas. Horizontal distance for the five- 
hour run is 90 kilometers. Bucket temperature (hourly) is solid 
line and radiometric (10.5 - 12.5 ym) temperature is dashed; 

„ other variables connect to their respective ordinate. The in- 

dicator isotherm is 22°C at 100 meters, which is approximately 
20 kilometers to the right of the maximum surface thermal gra- 
dient facing downstream. 
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to 1700 GMT is approximately 90 kilometers. 

Cross-section plots of surface salinity (hourly), chlorophyll-a 
(continuous; Lorenzen, 1966), volume scattering function, B 45 (hourly), 
radiometric surface temperature (dashed; continuous), and bucket temper- 
ature (solid; hourly) are drawn above the. cross-section of temperature 
down to 450 meters. Figure 2 is oriented such that the current flow is 
into the plane of the page and on the right-hand side (i.e., facing 
downstream). Surface observations are typical in that there is an in- 
crease of temperature and salinity, and a decreaseoof chlorophyll-a and 
scattering when entering the current from the cyclonic (left-hand edge 
facing downstream) side. The temperature field dominates the salinity 
to produce less dense water to the right in agreement with a geostrophi- 
cally balanced flow; chlorophyll-a and scattering interact to cause a 
color shift from greenish-blue outside the current, to the cobalt blue 
so characteristic of this relatively sterile water mass. 

Using the historical precedent of Leipper (1970) and others, the 
22°C isotherm at 100 meters depth was chosen as an indicator of the cur- 
rent. This isotherm is approximately 20 kilometers to the right of the 
cyclonic edge, and near the surface velocity maximum. The observational 
strategy greatly improved tracking efficiency as it gave the research 
vessel an average boost of 2 knots. The actual tracking was done using 
expendable bathythermographs (XBTs). After observing hydrographic sec-. 

tions across the Florida Straits from Cosgrove Shoal to Havana, and 

/ 

across the Yucatan Straits from Cabo San Antonio to Isla Contoy, a zig- 
zag tracking pattern was initiated. Typically, hourly BT's were taken. 


When the depth of the 22°C isotherm exceeded 125 meters or so, the ship's 
course was altered 90° to the left. This course was run until the 22°C 
isotherm was less than 80 meters or so, and then course was altered to 
the right. The pattern was continued from Yucatan, around the Loop to 
Dry Tortugas, in all but a few cases, where weather or fuel considera- 
tions made it advisable to run for Key West. After each trip, the posi- 
tion data were replotted and a smooth plot constructed, which made a 
best fit to the estimated courses and speeds made good. Positioning was 
accomplished using LoranA, radar, visual, and celestial observations. 

It is difficult to estimate errors, but based on Loran A-radar compari- 
sons, +1-3 kilometers seems reasonable. 

The expendable BT data were adjusted so that the average surface 
reading on the strip chart recorder was in agreement with the average 
bucket temperature. The strip chart recorder was tested (and adjusted if 
necessary) before and after each cruise with a test canister. Infrared 
temperatures were calibrated using a blackbody source in the laboratory 
so that the radiometric temperatures could be compared to the bulk 
^ values. This latter comparison is necessary in order to study the infra- 
red observations of the NOAA-2 scanning radiometer, and to properly 
interpret those data. 


TIME SERIES OF THE CURRENT 

d . 

Pathlines of the 22°C isotherm at 100 meters depth are summarized in 
figure 3. The series was obtained at 36-day intervals, coincidental with 
the 18-day period of ERTS. Dates of each survey are labeled on the 

} 
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Figure 3. Pathlines of the 22°C isotherm at 100 meters from August 1972, 
through September 1973. Surveys were 36 days apart and in 
coincidence with every other transit of the Earth Resources 
~ Technology Satellite: where two months are given, the survey 

extended over the change of date. Cruises were from three to 
six days in duration and the research ship motored with a 
following current. 
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appropriate pathline; where the cruise started at Yucatan in one month 
and ended off Dry Tortugas in another, both months are labeled. The 
shortest tracking time was three days and the longest six days, so that 
near synopticity was accomplished. Hydrographic station transects of 
the straits added two days to each cruise. Conditions prior to the 
first cruise were that a large anticyclonic detached eddy was observed 
in May, 1972. The current flowed essentially north from the Yucatan 
Straits and curved in a gentle anticyclonic curvature, ending tangent to 
the Florida Keys platform. 

By September conditions had changed markedly. The initial current 
direction had a significant easterly component and flowed directly toward 
the west Florida shelf. There was evidence of Loop Current water on the 
shelf, and the 22°C isotherm apparently went aground well north of Dry 
Tortugas. It is notable that a red tide of Gymnodinium breve was re- 
ported on the Florida shelf at this time. By late October the current 
had reformed to its southernmost extent, and evidence of Florida Bay 
water flowing south through the Keys was noted in both the ship track 
and an ERTS image. Murphy, et al_. (1974) have used this evidence to par- 
tially document the source of the first reported Florida east coast red 
tide. Apparently the organisms were carried by the current through the 
Florida Straits and exchanged with local coastal water north of Miami. 

By December 1972, the current had swung to the west and had pene- 
trated into the Gulf to the same latitude as August. At 24°N, the stream 
flowed in a sharp anticyclonic turn to the east. January 1973 was the 
only month in which transects of the Straits were not obtained. This 
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was due to 50-knot winds and high seas forcing the ship to turn back. 

Only four crossings were obtained but sufficient detail permitted 
the observation, for the first time, that the current penetrated north 
of Dry Tortugas (25°N). 

The "spring intrusion" (Leipper, 1970) continued through February, 
March, April -May, and June when the current extended to 27°N. As the 
current penetrated deeper it also swung further to the west and seemed 
to be controlled by the bottom topography of Campeche Bank in agreement 
with a potential vorticity conserving flow as postulated by Molinari and 
Cochrane (1972). North of 24°N the isobaths curve sharply to the west, 
and the current flows into deeper water where different dynamics probably 
dominate. . 

A deep intrusion in July 1973, coupled with a marked cyclonic curva- 
ture off the west Florida shelf, led us to expect separation of an eddy 
by the following cruise. The deep intrusion from the east was not a 
sampling artifact; a second vessel from the State University System of 
Florida Institute of Oceanography (SUSIO) obtained hydrographic station 
data across the shelf and out into the main current throughout this area. 
The furthest western extent of the current also occurred in July. 

By August 1973, the current system extended almost to the Mississip- 
pi Delta. The expected eddy had not separated. Very low salinity water 
(<24°/oo) was recorded by the SUSIO vessel, again cooperating with the 
NOAA ship,' all along the current-edge off the Florida Shelf. The hydro- 
graphic section in the Florida Straits had surface salinities less than 
30°/oo along the cyclonic boundary. Since the Loop Current was so close 
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to the delta area, and there are no other significant sources of fresh 
water, it seems probable that this was Mississippi River water in origin 

During the last cruise, September 1973, the current was found well 
to the south again, at approximately the same penetration as February. 

A trackline from Ft. Myers, west to 87°W, and north to Pascagoula, con- 
finned that an anti cyclonic eddy had indeed separated, and a significant 
change in the hydrography of the eastern Gulf had occurred in one month. 
There was no hint in the extensive August data that a recirculation had 
begun as a prelude to the eddy formation, although observations by 
Cochrane (personal communication, 1974) made between the April -May and 
June cruises, showed substantial closure in his isotherm field in this 
area. 

Niiler and Richardson (1973) obtained a time series of the trans- 
port of the Florida Current between Miami and Bimini. These data, 
taken during 1964-1970, show a transport minimum in December and a max- 
imum in June. The data in figure 3 show a minimum penetration into the 
Gulf in October-November and a maximum in July-August, which is in fair 
agreement with the average transport cycle. 

EDGE EFFECT 

The cyclonic edge of the Loop Current tends to concentrate flotsam 
and jetsam. Natural materials, such as surface marine algae, were seen 
from the June photography to have a pinkish cast on color IR (SO 397) 
film; this is confirmed in black and white IR photography (type 2424 fil 
and a 89B Wratten Filter). When these algae are present, . they appear 



as a bright lineation marking the edge of the current in ERTS MSS-6 
{0. 7-0. 8ym) . 

A second consequence of the boundary layer dynamics is to concen- 
trate chlorophyll bearing organisms. A typical profile across the cur- 
rent boundary is given in figure 4. This transect is from Key West 
harbor, out the channel into the coastal water, and across the front. 

The feature of interest is the peak in chlorophyll-a concentration just 
at the boundary. This is a phenomenon noticed frequently and occurs 
to a varying degree in the deep sea as well as near shore. Details of 
the variability will have to await further analysis (Ednoff, 1974) be- 
cause of expected seasonal and biological dependence. One would expect 
however. that a shift towards the green would occur and enhance the boun- 
dary in MSS-4(0.5-0.6um) . This does indeed happen as will be discussed 
below. 

A third observation is the usual change in sea state across the 
boundary. When winds and associated waves cross into zones of high cur- 
rent velocity with an opposing set, the seas generally build up rapidly. 
It is not uncommon for the sea state to increase from 1 to 2 meters when 
crossing into the current. Thus an increase of white caps, foam, and 
bubbles is encountered near the edge which increases the reflectance in 
the current. This raises the reflectance in all channels, but in a non- 
uniform {wavelength dependent) manner. 

The term "edge effect" (Maul, 1972) was coined to describe these 
phenomena which can be exploited to detect the boundary of. major ocean 
currents. Other streaming events, such as sediments being entrained by 
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Figure 4. Surface temperature trace (upper) and surface chlorophyll -a 
- _ profile (lower) across the Loop Current Front and into Key 

West Harbor, Horizontal scale across the figure is approxi 
mately 75 kilometers. 
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the Gulf Stream after passing source regions such as Cape Hatteras 
(c f. figure 13), can be considered part of this edge effect. It allows 
a multispectral approach to recognize the boundary of these currents in 
the absence of sea surface temperature changes. 

ERTS OBSERVATIONS OF OCEANIC FRONTS 

The location of the cyclonic boundary of the Loop Current was ob- 
tained by tracking the 22°C isotherm at 100 meters depth using expendable 
bathythermographs. The pathline of the isotherm is of the order of 20 
kilometers to the right (facing downstream) of the surface frontal zone 
(Hansen and Maul, 1970). This, coupled with ship observations of chloro- 
phyll-a, surface temperature, volume scattering function, and sea state, 
provides the baseline measurements to insure that the interpretations to 
follow are well founded. Figures 2 and 4 summarize the general conditions 
of increasing temperature and salinity and decreasing chlorophyll and 
scattering when crossing the boundary into the current and serve to 
orient them to the indicator isotherm. These observations confirm (Maul 
and Hansen, 1972) that changes in all these properties occur simultan- 
eously, allowing an increase in confidence when a recognition decision 
is made using several variables. 

Diffuse reflectance from beneath the ocean, which is defined as the 
ratio of upwelling to downwelling irradiance, measured just above the 
surface, is rarely more than 0.05. Reflectance from the ocean's surface, 
which is independent of diffuse reflectance, can be comparable or even 
substantially larger depending on sea state, and this has been shown by 



Maul (1973) to be a useful indicator of the current. Reflectance from 
clouds and agricultural scenes is sometimes an order of magnitude greater 
than from the ocean, even in the 500-600 nm wavelength (x) region. In 
order for the NASA Data Processing Facility (NDPF) to produce an image 
for an average scene radiance, the ocean signal is compressed into the 
lowest few gray scales. This is clearly illustrated in figure 5 which 
is a scanline plot across the boundary of the Loop Current from the 
mul ti spectral scanner. The large spikes in all four channels (MSS 4 
upper) are clouds; there seems at first glance to be very little change 
in digital number (DN), which is proportional to radiant intensity, as a 
function of the sample number. However, careful examination shows that 
the average value of the DN at samples greater than number 950 is slightly 
larger than those before this point. It will be seen that this marks the 
transition to higher radiances caused by increased sea state in the current. 


CONTRAST STRETCHING 

ERTS was not designed for ocean observations, however very early in 
the mission it became apparent that features of oceanographic significance 
were visible (e.g., Charnel 1 and Maul, 1973). In the coastal zone, ra- 
diance changes as large as 0.23 mW cm -2 sr“. 1 may be encountered across the 
interface of a waste disposal site, but in the deep sea, changes of less 
than 0.016 mW cm _2 sr _1 occur across a current boundary. Ross (1973) was 
able to demonstrate a degree of success with optical enhancement tech- 
niques using ERTS images in coastal areas, but changes whose radiance 
differences are more than an order of magnitude less- than those of waste 
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Figure 5. ERTS scanline plot across the Loop Current. front. Top scanline 
is MSS 4, next MSS 5, MSS 6 and MSS 7 on the bottom. The large 
energy spikes are clouds. At scan spot number 950 there is an 
increase in the average value of the digital number of 1 or 2; 
this marks the cyclonic edge of the current. 
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dumps cannot be optically enhanced because their signature is insuffi- 
ciently defined in the original film product. Any signatures sensed at 
all, however, will be in the digital data because ERTS digitizes the multi- 
spectral scanner (MSS) output onboard the spacecraft; such data are there- 
fore first generation information. Further, using the digital product 
permits actual radiances to be manipulated, thus avoiding certain vagaries 
associated with films, chemicals, and optical transforms. 

The full range of 128 digital number (DN) steps on MSS 4 (0.5-0. 6 urn) 
MSS 5 (0.6-0. 7 ym), and MSS 6 (0. 7-0.8 ym), and 64 steps on MSS 7 
(0. 8-1.1 ym) are required in order to digitize the range of radiances 
observed from the earth in each of these bandpass intervals. Figure 5 
summarizes the range of DN in the open ocean; other scanline plots show 
that DNs averaging 50 (-0.97 mW cm -2 sr -1 ) are typical over the land in 
MSS 4, but oceanic values in this channel average 20 (-0.39 mW cm -2 sr _1 ). 
In addition, the range of DNs over the land is much greater (ADN-25) as 
compared to the water (ADN~8). 

For the ocean scene then, all the information is constrained to a 
low level, and to a low range, of radiance. In order to optimize a dis- 
play of this information, a narrow range of DNs expanded about the mid- 
value of the ocean DN levels must be used. To do this a sub-matrix of 
pixels (called the training area), approximately 500 x 500 elements, of 
"cloud-free" ocean radiances are analyzed statistically for their mean 
(DN). standard deviation (o}, t and frequency distribution. If the training 
area is carefully chosen to be representative of the range of oceanic 
radiances, and the distribution is Gaussian, then 95% of all data are 
included in the range DN ± 2 $, 
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, To insure inclusion of all values in the original data, with em- 
phasis on a certain range, a mapping algorithm for the stretch variable 
( 5 ) must be chosen. If M is defined as the maximum value that the D/A 
converter can accept as an input voltage, a simple expression for con- 
trast stretching is: 



where n is a constant. The distinct advantage of a continuous function 
such as a hyperbolic tangent is that all values are included and their 
weights are adjusted according to the value of n. Several disadvantages 
such as nonlinearity of output and computational inefficiency favor a 
stretch variable of the form 


C = M, for DN < DN - 2a 


C = M 


3 ". 

9 9 


- 2o.<(DM - DN)< + 2a 


0 ) 


5 = 0 , for DN >_ DF + 2a. 

^ Equation 1 thus produces a negative, linear (for n = 1), statistically 
based contrast stretched image, of the data in the interval DN +_ 2o. A 
positive image could be produced if the numerator of the bracketed term 
in (2) were changed to DN - (DN - 2a), and the limits on z, = 0 and z, = M 
were reversed. 

Since ocean features are. 'always at the low end of the radiance range, 
it is sometimes advantageous to be able to emphasize them by a non-linear 
stretch. Equation 1 will allow this by setting n = 2, 3, etc. If the 
positive version of equation 1 is used, then n = 1/2, 1/3, etc. for 


stretching the "low DN values (at a great cost in computational efficiency) 
If a cloud-free training area is analyzed for its statistics, this objec- 
tive mapping algorithm allows a one-step computational effort. Using 
n = 1 is generally preferred since it produces a grayscale that is linear 
in radiance. In every case, the histogram from the training area's 
statistics should be studied to insure that the frequency distribution has 
a Gaussian form. If land or clouds were inadvertently included in the 
training area, or the area is too small, a non-Gaussian form will be 
encountered. 

Figure 6 is a negative print of an area due north of the Yucatan 
Straits using MSS 5 data; the scanline from figure 5 was taken through 
the center of this image. Computer enhancement in this image uses only 
five gray scales of the 128 levels available; all values below c = M are 
set to 0. The boundary between the resident .Gulf waters (left) and the 
current (right) is seen as a transition from light to dark tones respec- 
tively. Since the radiance levels in the ocean are so lovj the oceano- 
grapher-must resort to computer enhancement as this example shows. 

Figure 7 is another enhancement of the Loop Current boundary using 
MSS 5 where 7 <_ DN <_ 15 and n = 2. In this negative image of the western 
Florida Keys, water from Florida Bay extends into the Florida Straits and 
is entrained by the Loop Current. The current boundary in both figures 
4 and 5 was delineated by surface vessel tracks during the day of the 
satellite transit. Notice that the current is darker in tone (higher in 
radiance) in figure 5. This is caused by the dominance of -surface reflec- 
tiorf due to higher sea state in the current in figure 5, as compared to 
higher reflectance due to particles in the Florida Bay water in figure 7. 
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Figure 7. Negative print of computer enhanced (7 < DN ^ 15; n = 2) MSS 5 

image of Marquesa Key and Key West (ERTS ID 1099-15293). Change 
in radiance southwest of Marquesa from dark to light marks the 
ship located boundary between the higher intensity Florida Bay 
water and the lower intensity Gulf Stream. Bottom depth is in 
excess of 100 meters and thus does not contribute to the radiance. 
Horizontal distance across the image is 135 kilometers. 
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Figure 6. Negative print of computer enhanced ( 9 DN <^13; n = 1) MSS 5 
image of the cyclonic boundary of the Gulf Loop Current. Sur- 
face vessel track confirmed the location of the current to be 
the darker shade (higher radiance) region on the right hand 
_ side of the image (ERTS ID 1065-15411). Scanline plot in fig- 

ure 5 horizontally passes through the middle of the scene. 
Horizontal distance across the image is 90 kilometers. 
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For example, Ross and Cardone (1972) observed that with 15 m s~ 1 winds, 

10% of the surface of the sea is covered with white caps, and white caps 
reflect approximately 90%. If reflectance of the ocean is assumed to be 
5% in the absence of white caps and 10% of the sea is covered with white 
caps, then reflectance is R = 5% (.9) + 10% (.9) = 13.5% or 270% higher 
than without white caps. 

, SOME THEORETICAL CONSIDERATIONS 

At this point it is useful to consider in detail the processes con- 
tributing to the radiance spectrum N(^) at the position of the satellite. 
Solar radiation incident at the top of the atmosphere is absorbed and 
scattered in the atmosphere. Some radiation is scattered back into space 
without striking the ocean contributing a radiance N S U) at the satellite. 
The rest (which is not absorbed) will interact with the ocean. This inter- 
action can yield upwelling radiance above the ocean in three ways: 

1) specular reflection from the surface; 2) diffuse reflection from 
foam (bubbles) on or just beneath the ocean surface; and 3) the diffuse 
reflection from water molecules* and suspended particles in the water. Of 
these the first two phenomena are closely related in that they depend on 
the sea state, while the third source of radiance is essentially independent 
of sea state. The specular reflectance from the rough ocean surface can 
be computed by the methods of Cox and Munk (1954) in terms of the wind 
speed and incident radiance distribution. The diffuse reflectance from 
white caps can be approximately accounted for by assuming they are "white", 
Lambertian, and have an albedo of 1 so that they contribute uniform 

i 
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upwelling radiance just above the surface given by 



where f is the fraction of the scene covered by the white caps and H Q (\) 
the irradiance incident on the sea surface. The third source of radiance, 
that from beneath the surface, is the most difficult to compute and will 
be discussed in detail below. We can write the radiance at the satellite 
as 

N(X) = N S (X) + y(x)N ss (x) + ct(x) N d (x) (2) 

where N $s (x) is the contribution at the surface due to reflection from the 
surface and white caps, is the diffuse radiance just above the surface 
due to photons that have penetrated the ocean, : a(\) and y(x) are atmos- 
pheric transmittance factors for N d and N ss . a(x) and y(x) are in general 
not equal since the radiance distribution (variation with angle) of N ss 
and N d are different. It is possible for photons to be reflected from the 
surface, backscatter from the atmosphere into the ocean, and scatter back 
into the atmosphere. Photons which do this we consider to be a part of N d 
It should be stressed at this point that N d (x) in the above equation is 
the only source of radiance 'that contains information about conditions 
beneath the sea surface such as the concentration and composition of sus- 
pended particles and dissolved organic material. A thorough understanding 

r f 

of the dependence of N d U) on the basic optical properties of the water 
and. its constituents is required in order to obtain quantitative informa- 


tion about the constituents of the ocean from measurements of N d (x). 



The remainder of this section is devoted to relating N^{a) to these 
optical properties. We shall assume that the radiance distribution in- 
cident on the sea surface is given, and that the radiance is transmitted 
from the ocean to the satellite in a known manner given by a (x) and y (\) 
in equation 2. This reduces the problem to that of taking a known down- 
welling radiance distribution just above the sea surface, and computing 
the distribution (just above the surface) of the upwelling radiance. As 
mentioned above, radiant energy interacting with the ocean can be absorbed 
by water, suspended particles, and dissolved organic material commonly 
called yellow substance with absorption coefficients a w , ap, ay respec- 
tively and scattered by the water and particles with scattering coeffi- 
cients b w and b p (scattering by the yellow substance appears to be negli- 
gible). The total attentuation coefficient, c, for these interactions 
is given by 

✓ 

C = C„ + Cp + C y (3) 

where c w = a w + b w , 

c p = a p + b p , (4) 

c y - a y 

are the beam attenuation coefficients of the water, particles, and yellow 

substance. The scattering is further characterized by the phase function, 

P (e), which relates to the intensity of radiation, dJ(e), singly scattered 

d . 

from a small sample volume, dv, when illuminated by an incidence irradiance. 
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H Q , through 


P(e) = 


dj(e) 

(H 0 dv}b 



sine de = 1 . 


( 5 ) 


-The total phase function for water and particle scattering is 

p( 0 ) = (bwP w (e) + b p p p( e )) 

(b w + b p ) 


where P w and P b are the phase functions due to water only and particles 
only respectively. It is convenient to further define the forward (F) 
and backward (B) single scattering probabilities by 



and the single scattering albedo by 


Hence it is clear that 

B = (bpBp + b w B w ) 

<»p + b w> 

d (b w + bp) 

“o “ Ca w + a p "+'a y + b p + b w T * 
and in general , 0 <_ <_ 1 . 


( 7 ) 

( 8 ) 
( 9 ) 


f 
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It should be noted that aY\_ of the above quantities depend on wave- 
length (x). The transfer of radiation in the ocean is governed by the 
radiative transfer equation which has been discussed in detail by 
Chandrasekhar (1960) and Preisendorfer (1965), Gordon and Brown (1973) 
have computed the diffuse reflectance (upwelling irradiance/incident 
irradiance) just above a flat homogenous ocean as a function of its 
optical properties by a Monte Carlo technique. Using a combination of 
the parameters which arise naturally from the quasi -single scattering 
model (Gordon 1973), Gordon, Brown and Jacobs (in preparation) show that 
the diffuse reflectance (Rd ) can be written 

R d = 0.179X + 0.0510X 2 + 0.1710X 3 (10) 

Btoo 

where x = (V - „ 0 F) . 

R d is to first order independent of the distribution of the incident 
irradiance. It should be emphasized that the above equation does not 
include the irradiance specularly reflected from the sea surface or white 
caps, i.e. R d is the contribution to the reflectance from photons which 
penetrate the surface and are multiple scattered back into the atmosphere. 
Preliminary computations (Gordon and Brown unpublished) indicate equation 
9 is also valid for a moderately rough surface. The radiance distribution 
above the sea surface (due to photons scattered out of the ocean) is only 
weakly dependent on w 0 , and B'in the same fashion as R d . Hence, it is 
sufficient to study the influence of the optical properties of the ocean 
on R d alone. 
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Equation 9 shows that the important oceanic parameters are co Q and B, 
and so the observed reflectance spectrum, R^x), can be explained en- 
tirely through a knowledge of co 0 (x) or B(x). Conversely, under optimum 
conditions, only ^(x) or B(x) can be deduced from R<j and then only if 
one of these quantities is already known. The situation appears quite 
depressing when it is realized that to 0 ( x) is only imperfectly known even 
for pure water, and to our knowledge there are no measurements of Bp(x) 
and io 0 ( a) for various kinds of suspended particles. Hence, at the pre- 
sent time it is difficult to interpret to 0 (x) and B(x) even if both could 
be extracted from R<j(x) measurements. This underscores the necessity of 
laboratory experimentation to determine these optical properties for 
various ocean constituents such as marine phytoplankton, and suspended 
mineral particles, for adequate interpretation of oceanic "color". If 
the optical properties of the constituents are known, then it is theoret- 
ically possible to determine their concentrations through observations of 
Rd(x) as is discussed in several examples below. 

We shall now examine a particularly simple problem of interpretation, 
that of determining the concentration of suspended material in the absence 
and presence of yellow substance. The particles are assumed to be non- 
absorbing, and their scattering coefficient is assumed to be independent 
of wavelength [cases where bp = (const.) x~ n with n <_ 1 have also been 
investigated and yield results not dramatically different from the n = 0 

i 

case]. These calculations will not apply at all to locations with par- 
ticles containing absorbing pigments such as phytoplankton in the water, 
since Mueller (1973) and Gordon (1974) have shown that scattering from 



such particles varies strongly wi th wavelength near the pigment absorp- 
tion bands. 

In order to use equation 10, B(x) and w Q (x) in equations 8 and 9 
must be determined. We shall follow Tyler, Smith and Wilson (1972) and 
use for b w and B w values experimentally computed for sea water contin- 
uously filtered for 18 hours (Petzold 1972), rather than theoretical com- 
putations of these quantities for pure water. We do not feel comfortable 
with this, however, on the above basis they have estimated the absorption 
spectrum of "clear natural water", and since we use their a w (x) spectrum, 
this assumption is essential for consistency. It is thus assumed that 
B w = 0.1462 and b w = 0.008m -1 at 550 nm, so 



To find Bp, Petzold 1 s measurements of P ( 6 ) from the Tongue of the Ocean 
(T0T0), San Diego Harbor (SDH), and Off Shore -California (0SC) at 530 nm 
have been used. Assuming the above for B w and b w , the contribution from 
water is subtracted yielding the B p values listed in Table 1 (w 0 is given 
' in the parenthesis and w 0 for filtered sea water at this wavelength is 
0.136). 


TABLE 1 


Derived Values of Bp for the Indicated Locations 

' ■ .< ’ "" 


TOTO 

0.0165 (.59) 
0.0158 (.25) 
0.0130 (.26) 


SDH 

0.0186 (.82) 
0.0194 (.83) 
0.0169 (.91) 


OSC 

0.00966 (.59) 
0.00836 (.55) 
0.00836 (.55) 
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From Table 1 it is clear that B p not only varies considerably from one 
location to another, but also varies considerably for a single region. 
This is unfortunate in that for small w 0 , R d is directly proportional to 
(o 0 B, and w 0 and B are then equally important in determining R d for these 
cases. We shall assume that Table 1 gives the range of variation .of B p 
to be expected in natural waters, however this is probably not the case. 
Furthermore we will assume B p (and B w ) is independent of wavelength. 

Since a p is taken to be zero, and the a w (x) estimate of Tyler, Smith 
and Wilson 'is to be used, only a y (x) remains to be considered. This is 
taken from Jerlov (1968) and parameterized by 



where X is in nm. Reasonably high concentrations of yellow substances in 
the open ocean have ay = 5 x 10 -3 m _1 (i.e., near Galapagos Islands) how- 
ever, a v can be much larger in coastal regions. 

Using the above and equation 10, we have computed R d (x) called 
Reflectance (0,-) in the figures, as a function of b p /b w (proportional to 
the concentration of suspended particles) for various values of B p and ay 
The results are given in figures 8 through 12. 

Figure 8 shows the R d for x = 550 nm as a function of b p /b w and B p . 

The linearity of R d with b p /b w is to be noted. These results clearly 
demonstrate the importance in knowing B p for a quantitative determination 
of the particle concentration from R d (which is measured at the sea sur- 
face). Figure 9 shows spectra of Rd again for various values pf b p /b w but 
now with B d = 0.0165, These spectra have the same general shape as those 

r 

observed by Tyler and Smith (1970) just beneath the surface in Crater Lake, 
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Figure 8. Computed reflectance in percent at the sea surface as a function 
of the ratio of the particle scattering coefficient to the water 
. _ (only) scattering coefficient for 550 nm. The value of the 

fraction of backscattered light due to particles (Bp) for each 
curve is given on the right hand side. H 





A 

Figure 9. Computed reflectance in percent at the sea surface as a function 
of wavelength for various values of the ratio of the particle 
_ scattering coefficient to the water scattering coefficient. 

Mote the wavelength dependence in these spectra of changing 
the particle concentration only. 
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but are not in quantitative agreement even when the loss due to trans- 
mittance through the water surface is considered. Basically the dif- 
ference is that quantitative agreement requires too large a value of 
b p /b w . This could be due to the value-, of Bp used, the fact Bp is assumed 
to be independent of wavelength, or inaccuracies in their estimated a w (x) 
on which the present calculations are based. In any event the computa- 
tions can be used as a guide for examining ERTS related data. Figure 10 
shows the influence of yellow substance on the reflectance for the case 
with b p /b w = 128 and B p - 0.0165. The major influence of the Gelbstoff 
(yellow substance) is to depress the blue region of the reflectance 
spectrum. All of these results are summarized in figures 11 and 12, where 
the reflectances have been integrated over the ERTS MSS. channels 4 (fig-, 
ure 10) and 5 (figure 12) for various B p and a y as a function of b p /b w . 

It is seen that the Gelbstoff influences the reflectance in MSS 4 only 
at very high concentrations (for the open ocean) and essentially plays 
no role in the MSS 5 reflectance. Therefore, for. suspended particle 
concentrations, it appears that MSS 5 is best with MSS 4 only slightly 
' degraded by the dissolved organic materials. It should be emphasized again 
that the above discussion refers only to the case of no phytoplankton 
(chlorophyll)! 

We now turn to the problem of estimating bp/b w or some related quan- 
tity from satellite observations. Since the above calculations are only 
for the case of negligible phytoplankton (chlorophyll) concentrations, we 
must find criteria from which to choose which ERTS channels will satisfy 
this constraint. Also unless only MSS 5 is used we expect Gelbstoff to 
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REFLECTANCE (0,-)% 



Figure 10. Computed reflectance spectra, in percent, at the sea surface. 

The ratio of the particle to water scattering coefficients 
(bn/bw) and the fraction of backscattered light are 128 and 
0.0165 respectively. Values of the absorption coefficient due 
to yellow substance (ay) are listed on the left hand side. Note 
the shift of the spectral peak to longer wavelengths with in- 
creasing 

j 

. • ,/ 


- 32 - 




REFLECTANCE (O t -)% 



bp/b w 


Figure 11. Computed reflectance in percent at the sea surface as a function 
of the ratio of the' particle to the water scattering coefficients 
integrated over the spectral response of the MSS 4 filter (0.5 - 
0.6 ym). Values of the absorption coefficient due to yellow sub- 
" — stance (a y ) are listed along the right hand side. Note the lin- 

earity depends on both a v and the fraction of backscattered light, 

(B p ). 
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REFLECTANCE (0,-)% 



Figure 12. Computed reflectance at the sea surface, in percent, as a function 
of the particle to. water scattering coefficient ratio, integrated 
over the response function of the MSS 5 filter (0.6 - 0.7 ym). 
Value of the absorption coefficient due to yellow substance (ay) 

• - are listed along the right hand side. Note the linearity is not 

strongly dependent on ay but only on the fraction of backscattered 
light, Bp. 




be important especially in coastal regions where river runoff, etc., may 
be considerable. Furthermore there is the additional problem that the 
reflectance depends on Bp which is. also unknown, as well as bp/b w . We 
have developed a method which partially overcomes some of these problems. 
Considering the radiances observed in MSS 4 and 5, we have approximately 


N 4 “ “4 N d 4 + Y4N SS 4 + N $4 
N 5 = a 5 N d 5 + Y 5 N s$5 + .N $5 


(ID 


where a i and (i - 4, 5) are the fractions of N^- and N s$i (measured 
at the sea surface) that reach the sensor, It is assumed that a-j and yj 
are constant over an ERTS frame (if they vary in a known way, their in- 
fluence is easily accounted for and will not be discussed further). Now 
from the theory - 


N d 4 = k 4 5fi. 

b w 

N d 5 = k 5 b P 


( 12 ) 


where k^ and k^ are essentially independent of bp/b w , but depend directly 
on Bp. Taking the horizontal gradient, Vfj, of equation 11 using equation 
12 we find 

'hN, = TflV H N ss1 03) 


since N s . is nearly constant over a frame (vN*.. = 0). The last term in 

equation 13, v^N ss ^, is the horizontal gradient of the reflected radiance 

/ 
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from the sea surface. This is nearly zero everywhere except where 
the sea state changes dramatically with horizontal distance (for example 
in figure 6 ). Nearly everywhere on the frame then, v H N ss .. = 0, so 

v H N 4 . k 4 a 4 (14) 

v H N 5 k 5 a 5 

Again 134/05 is constant (or slowly varying) over the scene; so almost 
everywhere in the frame, variations in v^/v^Ng are the result of var- 
iations in k 4 /k 5 . Now if B is constant,, or if the "mean" B p is constant 
in each wavelength band, this would imply that essentially only the con- 
centration of scattering particles varies over the frame, i.e. not the 
nature of particles or their size distribution. In this case, k 4 /kg 
would be constant. Hence essentially k 4 /ks will vary, if the nature or 
size distribution of the particles varies over the frame (B changes), 

r 

or if the Gelbstoff concentration varies considerably over the frame, 
which would force to vary independently of kg. Thus if we find 


\ 


v H N 4 


5 H N 5 


= Const 


(15) 


it is reasonable to expect that only the particle concentration changes 
over the frame. In this case 

v h N 4 or v^Mg ~ (particle concentration) (16) 

d 

Note here. that these relations should also apply to scenes containing 
phytoplankton if they are the dominant scatterers; but for mixtures of 


f 
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phytoplankton and suspended white particles, one would expect equation 
15 to be violated over a scene if the relative concentrations vary 
drastically. To reiterate, if 15 holds, probably only the particle con- 
centration varies over the frame, and 16 can be used to measure the gra- 
dient- of its concentration. Since Bp is unknown, the actual concentra- 
tion is indeterminate without ground truth to better than a factor of 
two, since this is the assumed uncertainty in B p . Several examples of 
the use of the above ideas for analysis of ERTS data are presented below. 

DISCUSSION 

Equation 16 was applied to MSS 4 and MSS 5 data in figure 13. This 
computer enhanced negative image of the Cape Hatteras region shows what 
appears to be large gradients in suspended sediment. The Gulf Stream 
apparently has entrained particles from Rayleigh Bay and is carrying 
them out to sea. Along the scanline shown we have computed 

■§£ [N 4 ] = g|^- [ 1 N 5 ] + h (17) 

where the least squares values for g = 0.38 and h = 0.00, and the linear 
correlation coefficient, r,-is 0.66. The coefficient of determination, 
r 2 , which is the ratio of the explained variation from the mean (by the 
least squares line) to the total variation, is 44%. Physically this 
means that this image may not be useful for determining particle concen- 
trations because Bp, (i.e. ay), or bp/b w may be changing. The extent 
to which this holds true in natural waters is unknown; this will be 




Figure 13. Contrast stretched (4 < DN <_ 12; n = 2) negative image of the 

ocean area offshore of Cape Hatteras (ERTS ID 1132-15042). The 
Gulf Stream can be seen as the bright area to the south of the 
entrained sediment from the coastal estuaries. The least squares 
fit of eq. (17) was done along a scanline north of the Cape and 
extending from nearshore, through the suspended sediment and into 
the current. Extensions of this plume were observed for 150 
kilometers further east on other ERTS images. Horizontal dis- 
tance across the image is 135 kilometers. 
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extensively tested in the New York Bight area where turbidity measure- 
ments are being made by NOAA vessels concurrent with ERTS transits. A 
plot of the gradients from figure 13 and discussed above is given in 
figure 14. Note that MSS 4 follows MSS 5 quite well in this low 
passed filtered version of equation .16. 

The data of figures 9 and 10 can be compared to measurements of 
upwelling irradiance. Data given in figure 15 were observed using a 
1/4-meter Ebert spectroradiometer from 3 meters above the surface; these 
observations were made during the time frame of figure 7 and represent 
the water types shown in that image. All spectra were carefully selected 
to represent the same downwelling irradiance, sea state, sun angle, cloud 
cover, and absence of bottom influence. Specular reflection due to 
waves were minimized by preselecting ten spectra with similar shapes. 
After digitizing the records, averages and standard deviations (a) were 
computed at each wavelength; if values exceeded the average by la, they 
were rejected and a new mean computed. Absolute values of the spectra 
are traced to NBS through the 2-meter integrating sphere at NASA's 
Goddard Space Flight Center. 

As the chlorophyl 1-a concentration increases, the spectral peak 
shifts to longer wavelengths in similarity with the computed data in 
figure 10. This is not to say that the chlorophyl 1-a and the yellow 
substance produce quantitatively similar results; in general they will 
not. What it does imply is that the broad absorption bands in pigments 
such as chlorophyll and Gelbstoff at shorter wavelengths (<600 nm) 
produce qualitatively similar effects in the spectra. When other fac- 
tors are equal , increasing the concentration of pigments will cause a 
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decrease in the radiance in MSS 4 only. In nature however, such as in a 
plankton bloom or in a river plume, increased amounts of chlorophyll are 
normally accompanied by increased .particle concentrations (organisms 
which contain the chlorophyll) as well as yellow substances (decay pro- 
ducts of the plankton), and increased amounts of "white" suspended par- 
ticles in the river plume case. The radiances in MSS 4 and 5 will in 
these cases vary in a manner which will depend again on the relative 
concentration of the constituents, and since a and b vary almost inde- 
pendently, the reflectance signature, which depends on u> 0 (x) and B(x) is 
not unique and can be unraveled only when the variation of the optical 
properties of the individual components with wavelength is known. 

As an example of this, the spectra given in figure 15 were inte- 
grated over the MSS 4, MSS 5, and MSS 6 filter functions. A series of 
numerical tests were then made of ratios, differences, and sums to see 
if the three water types could be distinguished at the sea surface. It 
was quite easy to distinguish on the basis of such calculations between 
the Gulf Stream waters and the coastal waters, and between the coastal 
waters and the plankton bloom; however it was not possible to distinguish 
between the Gulf Stream and the plankton bloom. This suggests, as the 
theory implies, that for these data the rati o test (MSS 4/MSS 5) is not 
likely to be successful in specifying the chlorophyll-a concentration. 
Further, since the sea surface component, N $$ , (x), spectrally alters 
N(x), numerical tests (ratios', differences, etc.) are invalid indicators 
of oceanic properties. Theoretically, the signal of a differential 
radiometer (Arvesen, 1972) is subject to such variation in o) q . If it 
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Figure 15. Observed upwelling spectral irradiance in the Gulf of Mexico, 
November 1972. The three spectra represent typical observa- 
_ tions during the time series, and show the shift of the dom- 

inant wavelength to larger values with increased surface 
chlorophyll-a. The volume scattering function B (45 ) is for 
blue (436 nm) light. 
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is fortuitous that w always varies in the ocean such that R( 0 , - ) 

(J D / J 

is constant (Duntley, 1972), then a dual channel instrument will work 
in the absence of sea state changes. 

The solar spectrum is fairly flat in the visible region (H(0,-) 450 / 
H(0,-) 700 ~1 .3). Hence, higher sea states, which reflect high percen- 
tages of white light, proportionally add more long wavelength energy to 
the upwelling irradiance. In terms of a chromaticity diagram, this 
means that the purity changes but not the dominant wavelength. To cor- 
rect for sea state a channel in a multi spectral scanner near 1 pm appears 
useful. At this wavelength there is a maximum in the attenuation coeffi- 
cient of pure water (i.e. u o "^0), a maximum in atmospheric transmissivity, 
and the influence of chlorophyll seems low. 

Probably the most efficient method of determining the concentration 
of the constituents in the ocean will be to compare theoretical and ex- 
perimental spectra, adjusting the constituent concentrations in the 
theoretical spectra until agreement is found. This of course requires 
a basic understanding of the optical properties of the constituents which 
' ca n be derived only from careful in situ and laboratory experiments. It 
seems that at the present time much energy and money is being expended 
to try and use optical methods to locate and study materials with nearly 
unknown optical properties suspended or dissolved in a medium with only 
poorly known optical properties. This must be overcome before significant 
progress can be made. 

An interpretation of oceanic observations from ERTS is given in 
figure 16 which is a computer enhanced MSS 6 negative image of southeastern 
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Figure 16. Negative print of computer contrast stretched (7 £ DN 15; n = 
ERTS image of South Florida (ERTS ID 1026-15230). The Florida 
Current can be seen as a line of dark lineation parallel to the 
coast; bottom is essentially invisible in this MSS 6 scene. A 
ship can be seen by its characteristic V-shaped wake just off- 
shore of Miami Beach. Possibly the Virginia Key sewer outfall 
_ area can be observed by its low reflectance due to an organic 

slick. Horizontal distance across the image is 90 kilometers. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


- 44 - 




Florida. The dark lineation parallelling the coast in the upper portion 
of the image is a zone of high reflection caused by locally increased 
^ss a ^ on 9 ^e ec ^9 e of the Florida Current. The increase in surface 
reflectance is probably caused by surface wave interaction with the 
cyclonic boundary and is not bottom influence. This further explains 
the edge effect and supports the discussion on the local dependence of 
v rN on V H N SS only given above. Another example of the dominance of N s$ 
is the bright slick areas (low N) off the Virginia Key sewer treatment 
plant. This is probably -caused by the dampening of the glitter causing 
capillary waves in the oil film associated with the organic slick. The 
slick, which has drifted south past the popular Key Biscayne beaches, 
offers an explanation of: the narrow lineation off shore. in the Florida 
Straits: a passing oil tanker heading south which is pumping her bilges 

would cause a similar feature on the image. Thus ERTS could be useful 
in patrolling coastal waters for such illegal -acts which affect the near- 
shore water quality. 


NEARSHORE OCEANOGRAPHY 

The New York Bight is currently the subject of an intensive program 
of marine ecosystem analysis' (MESA) by NOAA. This program provides ex- 
tensive opportunity to relate basic oceanographic observations to the 
observed satellite radiances. Cruises by NOAA research vessels are 
being made in synchronization with ERTS transits. The transmittance 
data given in figure 10 were collected over a several day period nearly 
centered on the day of transit. 

t 
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Data can be smoothed before contouring in an effort to remove ex- 
cessive noise and high wavenumber variability. This is accomplished by 
convolving a boxcar smoother kernel (SK) with the data matrix (DM) 

SM = 3^IT ?? DM(j ’ k) {18) 

to produce a smoothed matrix (SM). 

Distinction should be made between smoothing and filtering. Smooth- 
ing in sample number space is defined as convolution with a boxcar filter 
kernel. Filtering in sample number space is convolution with a kernel 
which is the Fourier transform of a gate function from sample interval 
space. (Spectrally, sample number space is analogous to time, and sample 
interval space is analogous to period.) Smoothing is a form of filtering, 
but it employs a kernel that has undesirable side bands which can adverse- 
ly affect the output. 

In order to analyze the frequency characteristics of ERTS data, 
power spectra were computed along and normal to several scanlines in an 
oceanic region. Power spectra for data along a scanline (solid) and nor- 
mal to scanlines (dashed) are given in figure 17; 600 points and a Lanczos 
squared window were used to compute these results. 

Both power spectra have the shape typical of the geophysical con- 
tinuum, i.e., increasing power with increasing sample interval. However, 
the energy in the along scanline plot is generally higher than across the 
scanline. Strong (1973) argued that ERTS data have sunglint at the lar- 
ger nadir angles toward the solar specular region. That notion is sup- 
ported in these power spectra, since the energy at large sample intervals 
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is higher along the scanline {oriented toward the specular region) than 
across scanlines (which is quasi-parallel to the specular region). The 
large amplitude spike in the across scanline power spectrum of figure 17 
at every sixth sample is not observed in the along scanline power spec- 
trum. This spike is due to the so-called "six-line banding" caused by 
the use of six-detectors in the MSS. 

Many other power spectra were computed from different scenes and 
different spectral bands. There is a great deal of variability in the 
results because power spectra are very sensitive to the data chosen. 

Thus any conclusions based on Fourier estimates must be carefully viewed 
in light of the particular sample and judiciously applied to the ensemble. 

From study of representative power spectra such as these, a filter 
can be designed which will effectively eliminate undesired variations and 
noise. The Lanczos squared low pass filter used in this study is designed 
to eliminate the six-line banding. It has a half power point (6 db down) 
at 10 samples and is 90% effective (20 db down) at 14 and 6 sample points. 
The filter kernel was chosen to have 19 x 19 elements. That size is a 
compromise based on sharp cutoff for a 10-sample interval and a practical 
amount of computer time. 

Examples of different smoothing and filtering on the MSS 5 scene of 
the New York Bight are given in figure 18. The upper three panels are 
examples of using 3 x 3, 10 x,10, and 19 x 19 boxcar smoothers, top to 
bottom respectively. The lowest image is a filtered version using the 
19-x 19 low pass kernel derived from analysis of figure 17. After con- 
volving DM with the appropriate kernel, each image has been stretched 
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Figure 17. Solid line: Along scanline power spectra of oceanic data in 

MSS 5. Dashed line: Cross scanline power spectra of oceanic 

data. Power is in units of DN which has a radiance equivalent 
of mW cm“ 2 sr" 1 . 


Figure 18. Effects of smoothing and filtering on a portion of the data 

in figure 2 [b]. Upper [a]; 3x3 boxcar kernel. Second from 
top [b]: 10 x 10 boxcar kernel. Third from top [c] : 19 x 19 

boxcar kernel. Bottom [d]: 19 x 19 Fourier filter. The 

images are contrast stretched after convolution with the appro- 
priate kernel. Negative images. 
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for 7 DN 15 and n = 2 by use of equation 1. Smoothing clearly reduces 
the overall variability in the imagery, but the six-line banding, which 
has the appearance of corduroy material in the original data, is still 
present even with a 19 x 19 smoothing kernel. More significant however, 
is the loss of detail when smoothing as compared to filtering. Detail 
in the waste dump site (V shaped feature) for example, is retained in 
the filtered matrix (FM) as is variability in the plume region. The 
filtered version eliminates variability in the scene with wavelengths 
less than about one kilometer. 

Most of the contribution to filtering with a 19 x 19 kernel comes 
from the central 10 x 10 elements of the matrix. Comparing the 10 x 10 
SK image (18 [b] ) with the filtered version (13 [d]), note that varia- 
bility is retained in the plume region, and east of the waste dump site 
for the same stretched DN range in FM, but not in SM. 

Large sample internal variability, such as glitter, can be removed 
by high pass filtering. An alternate technique is to least squares fit 
a low order polynomial to the scanline and subtract out the contribution 
due to glitter. The least squares approach will not reduce the data 
matrix size due to convolution by a hiigh pass kernel, however, the phys- 
ical interpretation is more.-succinct when the data are viewed as being 
bandpassed in sample interval space. 

The images in figure 19 are computer enhanced, filtered, MSS 5 data 
in the apex region of the Bight using equation 1 and the 19 x 19 filter 
kernel described above. Transmittance is defined as the percent trans- 
mission measured by a 10 cm attenuation meter in the surface layer 




Figure 19. Percent transmittance for the apex of the New York Bight. 

Upper [a]: 29 August 1973 (ERTS ID 1402-15065) transit con- 
trast stretched and filtered as in 8 [d] . Lower [b]: 16 

September 1973 (ERTS ID 1420-15067) transit processed as in 
[a]. Oceanographic stations used are located by dots. Nega- 
tive images. Horizontal scale in each figure approximately 
50 kilometers. 
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(upper meter) of the water column. Isopleths of transmittance were con- 
toured from the ship data without regard to the satellite data for com- 
parison purposes. A cursory inspection of figure 19 [a] reveals that the 
lower transmittance areas are also areas of higher radiance (negative 
image has high radiance black),, but in 19 [b] this is not the case at 
all. The assignment of gray scales is arbitrary in each image. 

If N $s is constant (or high pass filtered out as required for the 
glitter problem) and if N s is also constant, then changes in N are due to 
changes in i.e,, changes in R d (0,-) as given by equation 10. In this 
equation, F, and hence B, depend on the particle concentration, the par- 
ticle size, as well as the size distribution of the particles. Due to 
differential particle settling, changes in all three variables relating 
to the phase function (p(e)/b) will occur in the Hudson River plume area. 
In addition, changes in a> Q can be expected here due to changes in the 
Gel bs toff distribution; these changes are independent of scattering and 
affect only absorption. 

In addition to the variables affecting the reflectance, non-synop- 
ticity must be considered. Semidiurnal bidirectional tidal currents 
dominate the flow in the entrance to lower New York Harbor, but offshore 
the current regime is rotary, superimposed on a highly variable but gen- 
erally net anticyclonic drift. Horizontal excursions for the rotary 
currents, are less than 5 kilometers, and net drifts are approximately 
6 cm sec -1 (-5 km day -1 ). Wi rids were light and variable during these two 
passes so there was no significant effect on the surface circulation. 

The transmittance data therefore represent a quasi steady state situation 
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especially in the region of the waste dump site and along the coast away 
from the harbor entrance. 

The inconsistency between observations of transmittance and ERTS 
radiances presented here is contrary to the conclusion by Klemas et al . , 
(1973) in that these data show both 'theoretically and observationally 
that agreement hinges on the behavior of u> Q and B. These experiments do 
support their conclusion on the usefulness of several analog enhancement 
techniques. Additional spectral studies. of R(0,-) and the ijn situ var- 
iables are now in progress in the Bight area. These measurements will 
give the proper background for interpreting ocean radiances from space. 
In light of these needs, filtered and contrast stretched images have the 
advantages of physical interpretability in terms of both the radiances 
and their spatial distribution, and are recommended for use in the ocean 
and coastal zone. 


SUMMARY AND CONCLUSIONS 

The spectral properties of the oceanic front associated with the 
Loop Current have been studied by ship and satellite observations and by 
radiative transfer theory. .It is seen that computer enhancement is re- 
quired to extract useful information from the ERTS data for the ocean 
scene. The current boundary can be detected by changes in the surface 

i 

reflectance, N ss , as well as ‘the diffuse reflectance, , from below the 
surface; Nj however is dependent on both B and w 0 , and thus the spectral 
interpretation of ocean color requi res surface truth measurements for 
meaningful results. Particle concentrations, which can delineate currents, 
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can be estimated in MSS 5 if the ratio v h n 4 /v H% is reasonably constant 
over a scene. MSS 4 is strongly influenced by yellow substance, and par- 
ticle estimation based on these data are invalid in many coastal zones. 
Water mass identification using ratios or differences of MSS 4, 5, and 6 
data have no validity in either theory or observation in the Gulf of 
Mexico. Finally, it must be emphasized that the spectrum of upwelling 
radiance just above the surface is a function of both N d and N $s and that 
N s$ frequently dominates. 
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APPENDIX 


In this section several ERTS-1 images of the Gulf Stream System are 
summarized. Each of the selected images used in this report are located 
in. figure A.l. Frame numbers 1, 2, 3, and 5. were discussed in the body 
of the text and will not be repeated here. The selected images are from 
the area of concern as stated in the contract and does not mean that 
these were the only good images of the current. About 10% of the images 
received, that passed the 60% cloud cut-off criterion, actually had ocean 
current information. The rest were either in the wrong area or the cur- 
rent was not visible; the latter cause accounted for about 25% of the null 
information. 

Figure A. 2 was used by Maul and Gordon (1974) to study the capability 
of ERTS to determine the particle concentration. In that study, they 
showed that if 


^!cc c “ constant 
vMSS 5 


A.l 


then 


VMSS 4 or 5- vb p /b w A. 2 

which is in turn proportional to the gradient of the particle concentra- 
tion. Along the scanline shown in figure A. 2, the gradient in MSS 4 and 
MSS 5 has been computed using several filters; the correlation coefficient 
(r) has been computed for each filter pair as well. The results of using 
different low pass filters and filter lengths are summarized in figure A. 3 



The general trend is that r increases with longer filter kernels and 
with lower pass filters. This means that the gradient on the earth's sur- 
face with 8 km (100 sample low pass) spatial scale, has greater coherence 
between bands. It is not certain as to why the small scale correlation 
coefficient is so low. This area has significant amounts of coastal sedi- 
ments which are rich in organics, and thus the Gelbstoff concentrations 
probably vary over a wide range. 

In the images to follow, the interpretations given are based on the 
discussion in the text. Only one band is given in the appendix, but all 
were studied for interpretation. Clouds were identified in MSS 6 and 
MSS 7 and were optically masked out of the image before study, so that the 
data discussed are water phenomenon and not an atmospheric contribution. 
Water color shift from green to blue was made by studing if the variation 
in MSS 4 disappeared in MSS 5; if so, the color of low radiance scanspots 
in MSS 4 must be bluer than the higher radiance data. Similarly sediments 
are not seen in 'MSS 7 or frequently even in MSS 6; clouds are always seen 
in all four channels. The images are presented for the general use of 
other marine scientists without .discussion except for the figure captions. 




Figure A.l Location map of selected images in appendix A. Images 1, 2 , 

. ^ 3, and 5 were discussed in the text and are not repeated here. 

Images 4, 6-17 are labeled figures A. 4, A.6-A.17; there is no 
A. 5 and A.1-A.3 are used for introductory discussion. 
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Figure A. 2 Contrast stretched (4 _< DN 12; n=2) negative MSS 5 image of 
the ocean area offshore of Cape Hatteras (ERTS ID 1132-15042) 
observed on 2 December 1972. The Gulf Stream can be seen as 
the bright area to ,the south of the entrained sediment from 
the coastal estuaries. The least squares fit of eq. (17) was 
done along the scanline north of the Cape and extending from 
_ nearshore, through the suspended sediment and into the current. 

Extens-ions of this plume were observed for 150 kilometers 
further east on other ERTS images. Horizontal distance across 
the i.mage is 135 kilometers. 





Figure A. 3 Gradient of radiance received at the satellite in MSS 4 along 
the ocean line shown on figure A. 2. Scan spot numbers arbe- 
trarily start at 1 which is just off the image in the previous 
figure and over the coast. Correlation coefficient (r) for 
the MSS 4-MSS 5 linear correlation is: upper - 10 sample low 

pass (199 element kernel), r=0.085; middle - 100 sample low 
pass (59 element kernel), r=0.398; lower - 100 sample low pass 
(159 element kernel), r=0.698. Compare with figure 14 for 100 
sample low pass using 199 elements in the filter kernel, where 
r=0. 740. 
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Fiqure A. 4 Contrast stretched (16 <_ DN 5 24; n=l) negative MSS 4 image of 
eddies in the eastern Gulf of Mexico (ERTS ID #-1153-15292) 
observed on 23 December 1972. The eddies were visible in only 
one band and therefore are minus green in color. This implies 
that they are spin-off eddies from the Gulf Loop Current that 
have come onto the west Florida Platform where they will even- 
tually spin down exchanging momentum, nutrients, and heat with 
that water mass. Horizontal distance across the image is 45 
km which means the eddies are 20 km. Sense of rotation is un- 
known but probably anti cyclonic. 
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Figure A. 6 Contrast stretched (12 <_ DN < 21; n=l) negative MSS 5 image of 
Miami, Florida (ERTS ID 1350-15231) observed on 8 July 1973. 
Government Cut arid the Miami harbor basin were being dredged 
when this image on the ebb tide was taken. The dredge spoils 
can be seen as a jet flowing out into the shelf where a cyclon- 
ic circulation is caused by the Longshore drift. A cyclonic 
eddy, either spin-off or shear induced (Lee, 1974) by the Gulf 
Stream is seen in the upper right due to changes in surface 
reflectance. The Gulf Stream front at the lower right is seen 
due to higher radiance in the current. Complications of near- 
shore currents can be deduced because the offshore eddy and 
the northern inlets have a southward inshore drift that has 
reversed or slackened near the dredging plume. Horizontal 
distance across the image is 45 km. 
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Figure A 7 Contrast stretched sequence of MSS 6 data off Miami Beach, 
Florida (ERTS ID-1026-15230) observed on 18 August 1972. 
Sequence shows the effect of varying the DN range and n in 
the contrast equation, a) Upper left: Raw data, positive 

print, b) Upper middle: 7 < DN <J5; n=l . c) Upper right: 

7 < DN < 15;' n=3. d) Lower left: 7 <. DN <. 11 ; n=l . e) Lower 
middle: 9 < DN £ 13. f) Lower right: MSS 6 minus MSS 7 

then 5 ^ DN £ 13; n=3. Note that using powers of n brings 
out water details (cf. b & c) as does lowering the upper 
limit on DN (cf. b & d). Water details in e are less distinct 
than in b which uses the statistically determined limits. Sea 
state only is observable in MSS 7 so MSS 6-MSS 7 should mini- 
mize the surface contribution as is shown in f. Horizontal 
distance across each image is 45 km. 
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Figure A. 8 Contrast stretched (9 _< DN £ 17; n=l ) negative MSS 5 image of 
the region offshore of Cape Canaveral, Florida (ERTS ID 1260- 
15233) observed on 9 April 1973. Fine scale "shingle structure" 
(von Arx, Bumpus, and Richardson, 1955) can be seen along the 
cyclonic edge of the Gulf Stream. The shingle structure in 
this image is at least an order of magnitude less than von Arx 
et al_. description off the Carolinas. Wind streaks and a ship 
wake can be observed in the Gulf Stream which is leaving the 
coast at the lower right. Loss of water detail in Lake 
Okeechobee is caused by the contrast stretch limits and shows 
_ that the very high radiance from this water body is due to a 

turbidity about equal to the coastal inlet plume-in the upper 
left. Horizontal distance across the image is 135 km. 
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Figure A. 9 Contrast stretched (8 <_ DN £ 16; n=l) negative MSS 5 image of 
the Florida coast south of Cape Canaveral (ERTS ID 1206-15232) 
observed on 18 August 1972. Plumes of water from behind the 
barrier islands can be seen ebbing through the inlets. Com- 
plicated patterns of coastal sedimentation show a multifronted 
boundary between the nearshore water and the Gulf Stream. 

- Changes in the radiance appear to be dominated by changes in 

suspended sediments. Horizontal distance across the image is 
45 km. 
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Figure A. 10 Contrast stretched (7 _< DN j< 15; n=2) negative MSS 5 image of 
Cape Lookout, North Carolina, (ERTS ID-11 15-15152) observed 
on 14 November 1972. The shoals offshore of the promitory 
are shallow and keep sediments in suspension. By comparing 
the sediment pattern over such areas (cf. A.ll) an estimate 
of the bottom topography can be made, and conversely, the 
variability of the suspended load as well. There is no 
evidence of the Gulf Stream front in this image. Horizontal 
distance across the image is 135 km. 
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Figure A. 11 


contrast stretched (4 < DN < 1 1 • n-i \ 

Onslow Bay, North Carolina Terts jn i ?ns 9 i ci'rox MSS 5 image of 
13 February 1973. ComDlex n*tl» ID ~ 1 105- 1 5] 53 ) observed"™ 
to be well correlated with the "bottom** turtnd ity do not appear 
Pixels near the coast are an artifa?? ^ P ?F raphy ' The white 
routine wherein the NOAA computed ass inL** 6 contrast stretch 
the upper cutoff value (DN=11 in th?! 9 x ne 9ative zero to 
zero is reformatted onto h d oih f C3Se) ; K ' hen the negative 
vertently becomes a maximum v.nl I- u“ analo 9 ta P e » it inad- 
as a white pixel. The s?eL Whlch appears on the image 
and vvhite (DN=11) values are infp • 9 ? y because black (DN>11) 
^ images with n>2, h s b IT ^ , This do « not appear 
scale distance across the image is“kT£? ^ 


y 
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Figure A. 12 Contrast stretched (5 < DN < 13; n=2) negative MSS 5 image of 
the Cape Hatteras, North Carolina area (ERTS ID-11 32-1 5094) 
observed on 1 December 1972. Mairs (1970) correlated the 
pulses in the Ocracoke Inlet sediment plume to several stages 
in the tide. The Cape Hatteras materials are transported 
from further north, and an eddy has been observed in several 
ERTS images in the Hatteras region. Compare again the sedi- 
— ment suspension over Lookout Shoals with the scene in figures 

A. 10 and A. 11; detail changes are due to both different stretch 
limits as well as environmental conditions. Horizontal dis- 
tance across the image is 90 km. 
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Figure A. 13 Contrast stretched (14 £ DN £ 22; n=2) negative MSS 4 image of 
an eddy over Diamond Shoals, North Carolina (ERTS ID 1186- 
15093) observed on 25 January 1973. A deep penetration of 
water from north of Cape Katteras can be seen as part of the 
cyclonic eddy developing over the shoals. The Gulf Stream 
front appears to be well offshore of the eddy front. Meander 
scale in the Gulf. Stream appears to be about 50 km with an 
amplitude of about 10 km. Details in the sediment suspension 
patterns supports the notion of multiple fronts, probably 
tidal ly pumped and advected by alongshore current patterns. 
Streaks can be seen in the current at the lower right and may 
be internal wave activity. Horizontal distance across the 
image is 90 km. 
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Figure A. 14 Contrast stretched (7 < DN < 15; n=2) negative MSS 5 image of 
a sediment filament in the offing of Cape Hatteras, North 
Carolina (ERTS ID 1222-15093) observed on 2 March 1973. The 
squaring in the stretch equation over-accenuates the low 
radiance features, so that details at Cape Hatteras are lost; 
generally n=l is a better choice because it is linear with 
radiance for spectral comparison purposes. The sediment 
filament extends 120 km offshore and is probably material 
from Pamlico Sound as Ford, Longard, and Banks (1952) have 
suggested in their entrainment studies. Horizontal distance 
across the image is 135 km. 
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Figure A. 15 Contrast stretched (16 £ DN £ 24; n=3) negative MSS 4 image 
of the Gulf Stream in deep water off Cape Hatteras (ERTS ID 
1221-15035) observed on 1 March 1973. Taken one day earlier 
than A. 14 this plume is over 180 km long on this image and 
extends over 300 km from Cape Hatteras. The cumulus band is 
a frequent feature of the stream from space and is caused by 
the convection when cool dry continental air is heated by the 
stream. This feature follows the Gulf Stream front quite well 
as seen in the image and can be used as an ancillary feature 
for identification of the Stream. Horizontal distance across 
the image is 90 km. 
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Figure A. 16 Contrast stretched (12 <_ DN <_ 1*6; n=l ) negative MSS 4 image of 
a sediment filament offshore of Cape Hatteras, North Carolina 
(ERTS ID 1149-15033) observed on 19 December 1972. This 
image was observed 17 days later than the data in A. 2, but is 
an almost perfect extension of that sediment filament and ex- 
» tends about 270 km to sea. Note the DN limits have a range 

of only 5. The Gulf Stream waters are at the same intensity 
„ as the waters on the north of the sediment filament and the 

Stream is observable only because little small scale mixing 
has occurred across the boundary. Horizontal distance across 
the image is 135 km. 
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